Anomalous Flux Flow Resistivity in Two Gap Superconductor MgB 2 
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The flux flow resistivity pf associated with purely viscous motion of vortices in high-quality 
MgB2 was measured by microwave surface impedance. Flux flow resistivity exhibits unusual field 
dependence with strong enhancement at low field, which is markedly different to conventional s-wave 
superconductors. A crossover field which separates two distinct flux flow regimes having different 
Pf slopes was clearly observed in -BT|| ab- plane. The unusual //-dependence indicates that two very 
differently sized superconducting gaps in MgE>2 manifest in the vortex dynamics and almost equally 
contribute to energy dissipation. The carrier scattering rate in two different bands is also discussed 
with the present results, compared to heat capacity and thermal conductivity results. 



PACS numbers: 74.25.Fy, 74.25. Nf,74.70.Ad 



Since the recent discovery of high temperature super- 
conductivity in MgB2, great attention has been directed 
towards understanding the detailed nature of supercon- 
ductivity Q, Q- MgB 2 consists of two bands having 
roughly equal density of states (DOS), strongly two di- 
mensional er-band (sp x p y ) and three dimensional 7r-band 
(p z ) bands. Associated with these two bands, two very 
differently sized superconducting gaps are quite distinct 
and manifest themselves clearly in many physical prop- 
erties, which have become one of the central topic of 
the physical properties of MgB2- In fact, apart from 
the much debated pairing mechanism of high transi- 
tion temperature superconductivity, physical properties 
in the superconducting state i.e., the distinguishing char- 
acteristics of a superconductor with two gaps, are still 
largely unexplored. Recent scanning tunneling micro- 
scope (STM) 0, heat capacity Q, and thermal conduc- 
tivity [5j measurements have revealed that the physical 
properties in the vortex state of MgB2 are dramatically 
different from those expected in conventional s-wave su- 
perconductors. Despite these extensive studies, however, 
the dynamical properties of the vortices, such as the en- 
ergy dissipation associated with the vortex motion which 
is intimately related to the electronic structure [6|,lj]>l2|>l2|> 
are still far from being completely understood. 

When a vortex line in a type-II superconductor re- 
sponds to a driving current, frictional force is given by 
the damping viscosity, which in turn depends on the en- 
ergy dissipation process of quasiparticles in and around 
the vortex cores. To gain an understanding of this en- 
ergy dissipation, experimental determination of the free 
flux flow (FFF) resistivity is particularly important. FFF 
refers to a purely viscous vortex motion of vortices, which 
is realized when the vortices move in pinned free states. 
The FFF state in conventional s-wave superconductors 
has been extensively studied and a rather good under- 
standing has been developed. In s-wave superconductors, 
quasiparticles trapped within the vortex core that form 



Caroli-de Gennes-Matricon bound states play a key role 
in the dissipation process 0,0- The FFF resistivity pf 
is essentially proportional to normal state resistivity p n 
and applied field, 
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This Bardeen-Stephen relationship has been confirmed in 
most diry and clean s-wave systems throughout almost 
the whole Abrikosov phase. Recently, striking deviation 
of FFF resistivity from Eq.(l) has been reported in un- 
conventional superconductors with gap nodes 0, EH 
llll Il2l Il3|. For instance, in d-wave Bi2Sr2CuOg+5, Pf 
exhibits a //-linear dependence at low fields, followed by 
a non-linear //-dependence described well by, 
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at higher fields [9j- This demonstrates that the energy 
dissipation is strongly influenced by the superconducting 
gap structure. However, despite these extensive studies 
of the FFF state, the detailed microscopic mechanism of 
energy dissipation is not well understood, exposing our 
incomplete knowledge of vortex dynamics. This situation 
thus calls for a new textbook example of FFF resistivity 
in a superconductor with a distinctive gap structure. The 
FFF resistivity of MgB2 with two distinct superconduct- 
ing gaps may provide an unique opportunity for studying 
how the energy dissipation occurs when vortices move in 
the superfluid. 

In this paper, FFF resistivities of MgB2 are exam- 
ined by measuring microwave surface impedance. FFF 
resistivity exhibits unusual field dependence, which is 
markedly different from conventional s-wave supercon- 
ductors. These results are discussed in terms of energy 
dissipation in the superconductors with two distinct gaps. 

High quality single crystals of MgB 2 (7^=38.6 K) were 
grown using a high pressure method |14j . Typical size 
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was 0.3 x 0.3 x 0.1 mm 3 . Microwave surface impedance 
Z s = R s + iX s , where R s and X s are the surface re- 
sistance and surface reactance, respectively, was mea- 
sured by the standard cavity perturbation technique |15| . 
Cylindrical cavity resonators made of Cu (Q ~44000) 
and superconducting Pb (Q ~ 1.0 x 10 6 ) operated at 
28.5 GHz were used in TEon mode. Sample was placed 
in an antinode of the oscillatory magnetic field H u . Z s 
was measured in two different configurations, H u \\H\\ c- 
axis (Z^ 1 "" ) and \\H\\ a6-plane (Zj ? "'" ab ). 

Figures 1 (a) and (b) show the T-dependence of in- 
plane and out-of-plane surface impedance, Z^ b — R^ b + 
iX® and Z c s = R c s + iX^, respectively, measured by 
Pb-resonator in the Meissner phase. Both sets of data 
were obtained by the measurements of zf" " c and Zf" " afc 
(see the insets of Fig. 1). In the configuration H u | c, 

Zg b = Z^" c because the oscillatory currents flow within 



the ab- plane. For H u \\ ab, Z, 



H„\\ab 



is determined by 

H u \\ab _ 



the geometrical mean value of and Z% as Z, 

S g£b+gc Zs ■ Here S ab and S c are the areas of the faces 
where the screening currents flow in the a6-plane and 
along the c direction, respectively 0. Since S ab ~ 35 c , 
and Z c s ~ X.lZf' from the anisotropy of normal state 
resistivities (p c J p ab ~ 3 [13), Z c s and Zf make contri- 
bution in the same order to z^^ ab . In both configura- 
tions, both R s and X s decrease rapidly with decreasing 
T below T c . In the Meissner phase, microwave response 
is purely reactive and R s ~ and X s = /iowA, where po 
is the permeability of free space, lu/2it is the microwave 
frequency, and A is the London penetration length. On 
the other hand, response is dissipative in the normal state 
and R s = X s = p Q u>5 n /2, where S n = ^/2p n / p^w is the 
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FIG. 1: (a) Temperature dependence of in-plane surface 
impedance, Zf = Rf + iXf b , in MgB 2 single crystal. (b)The 
same data for out-of-plane surface impedance, Z C B = Rg+iX^. 
pt = 2pQ cm, p c n /pl b = 3 and S ab = 3S C were used. 
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FIG. 2: (a) _ff-dependence of surface resistance i?;f"" c nor- 
malized by the normal state values measured in H\\ || c 
at 40 K (diamonds), 30 K (inverted triangles), 20 K (trian- 
gles), 10 K (squares), and 2.3 K (circles). Inset: A plot of 
impedance R s /X s {0T) vs X s /X s (0T) for H\\ \\ c (filled 
circles) and for H\\ \\ ab (open circles). The solid lines 
are the results calculated from Eq. (3) at different pinning fre- 



measured in H II H 



quencies. (b) The same plot for R, 

ui || ab. Inset: T-dependence of i/° 2 ° and i?c2, as determined 
by the field at which R s is the normal state value. 
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skin depth and p n is the normal state resistivity. The 
absolute value of R\ and X l s (i = ab and c) were de- 
termined by comparison with dc resistivity p l n assuming 
that R\ = X\ = \J poUjp l n /2 in the normal state. Using 
p°^ ~ 2p£t cm and p^j p°^ ~3 at the onset, the in-plane 
and out-of-plane penetration lengths, \ a b(T — 0) ~ 500A 
and A C (T = 0) ~ 750 A, respectively, were determined. 
Our \ a b is smaller than the reported values jl8j . but 
the iJ-dependence of Z s is little affected by the abso- 
lute value of X a b- The obtained value of X c /X a b ~ 1-5 is 
much smaller than the anisotropy of upper critical fields, 
H ci/ H c2 - 4 ( see the inset of Fig. 2(b)), where H^ b and 
H^ 2 are the upper critical field in H parallel to the ab- 
plane and c-axis, respectively. We note that the small 
anisotropy of the penetration lengths is consistent with 
the recent result in Ref.|l9|. 

Figures 2 (a) and (b) show independence of i?^"" c 
( H u \\H\\ c) and R^ ab ( H u \\H\\ ab) measured by 
Cu-resonator. In both configurations, R s increase rapidly 
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with H . In the vortex state, Z s is determined by vortex 
dynamics. According to Coffey and Clem, Z s can be 
expressed as , 



Z s = i/iQiuX 
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where 5?. 
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is the pinning fre- 



quency, Sf = y/2pf/[iQLu is the FFF skin depth, S qp — 
y / 2p qp / pquj with the quasi-particle resistivity p qp is the 
normal- fluid skin depth |20j |. To estimate the pinning 
frequency of MgB 2 , plots of X s vs R s in two configura- 
tions are shown in the inset of Fig. 2(a). For comparison, 
the same data at several uj/uj p calculated from Eq.Q are 
also plotted. The inset of Fig. 2(a) strongly suggests that 
Lu/ujp is much larger than unity in both configurations, 
demonstrating that the FFF state is nearly realized in 
the present frequency range. In the FFF state, two char- 
acteristic length scales, namely A and Sf, emerge in ac- 
cordance with microwave penetration. In low fields, A 
greatly exceeds 5f(X 3> Sf). In this region, R s and X s 
are given by R s ~ P//A and X s ~ p.QU>\. In high fields 
where St exceeds A(<5/ 3> A), the response is similar to 
the normal state (R s ~ X s ), except that S n is replaced 
by Sf. In MgB2, this cross over occurs at a very low 
field (p,oH < 100 mT) in both configurations. This in- 
dicates that the influence of penetration lengths on the 
vortex dynamics is negligibly small except at very low 
field regimes. 

Having established that the FFF state is realized in 
the present experiments, FFF resistivities in H paral- 
lel to the c-axis and a6-plane, p^" c and p^" a& , respec- 
tively, are estimated next from Figs. 2(a) and (b), assum- 
ing uj/lu p = co. Figure 3 (a) shows pf and pf^ ah , 
normalized by the normal state values as a function of 
H/H c2 at T=2.3 K (T/T c =0.06). For comparison, the in- 
dependence of the FFF resistivity of single gap s-wave su- 
perconductor, as given by Eq. QJ, and of d-wave Bi:2201, 
as reported in Ref. 9| are also shown. The field depen- 
dence of both pf^ c and pf^ ab are spectacularly different 
from those of s-wave and d-wave superconductors. Both 
pf^ c and pf^ ab increase steeply with H at low fields, 
even steeper than pf in d-wave superconductor. In par- 
ticular, at H/H^ ^ 0.1, pf l|ab of MgB 2 is nearly half of 

pf^ ab at H C 2, indicating a striking enhancement of en- 
ergy dissipation at very low fields. This unusual FFF 
state of MgB2 can be seen more clearly in Fig. 3(b), in 
which the same data is plotted as a function of ^H/H c2 . 
A crossover field H cr which separates two distinct flux 
flow regime with different pf^ ab slope is distinguishable 

at y H cr /H^2 — 0.3, as indicated by the arrow. On the 
other hand, such a crossover cannot be easily identified 

H\\c 

mp f . 

The appearance of the crossover field can be attributed 
to the two distinct superconducting gaps. Unusual elec- 
tronic structure in the vortex state of MgB2 has been 
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FIG. 3: (a) Flux flow resistivities pj, (solid triangles) and 
pH\\ab ^ so jj ( j c i r cles) at 2.3 K normalized by the normal state 
values as a function of H normalized by the upper critical 
fields. The solid line represents the FFF resistivity of con- 
ventional s-wave superconductor expressed by Eq. (1). The 
dashed line represents the FFF resistivity in d-wave super- 
conductor reported in Ref.|9j. (b) The same data plotted as 



a function of v/ H/H c 



reported in several experiments. According to STM mea- 
surements, the vortex radius is much larger than ex- 
pected from H C 2 and superconductivity is strongly sup- 
pressed by a low field. Sharp increases in the heat capac- 
ity Ce H and thermal conductivity n e [f| at low fields 
with initial slopes much larger than in conventional s- 
wave superconductors also revealed the unusual vortex 
state of MgB2. In particular, in H || ab, C e is nearly half 
of the normal state value C™ at H/H^ ~ 0.1, followed 
by a gradual increase up to H C 2- Moreover, n e in H\\ ab 
becomes nearly ii-independent after n e reaches nearly 
half of the normal state value n n at H/H^ ~ 0.1. These 
strong enhancements of C e and n e at low fields has been 
attributed to the suppression of the small superconduct- 
ing gap originating from the three dimensional 7r-band. 
Therefore, it is natural to consider that the steep increase 
and crossover behavior observed in the FFF resistivities 
could also be attributed to the two superconducting gaps. 

In conventional s-wave superconductors, the energy 
dissipation associated with FFF is caused by the quasi- 
particles trapped within vortex cores 0] . Since in con- 
ventional s-wave superconductors the DOS within cores 
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is approximately ii-independent, FFF resistivity is pro- 
portional to the number of vortices as given by Eq.Q). 
The situation in MgB2 is very different. Below H cr , 
quasiparticles are trapped by both small and large gaps. 
Since the coherence lengths associated with the small gap 
are much larger than those associated with the large gap, 
the number of quasiparticles trapped within cores well 
below H cr is much larger than the number trapped solely 
by the large gap. This gives rise to strong enhancement 
of FFF resistivity at low fields. As H approaches H cr , the 
small gap is strongly suppressed. As a result, the contri- 
bution of bound quasiparticles of the 7r-band essentially 
becomes negligible above H cr . This gives rise to a gen- 
tler increase in FFF resistivity above H cr up to H C 2 . The 
kink structure observed at H cr in Fig. 3(b) implies that 
this crossover occurs in a narrow field range. Similar to 
the d-wave superconductor described by Eq.(J2J), pj^ ab 

below and above H cr exhibit nearly y/~H dependence. To 
clarify the origin of this if-dependence, more detailed in- 
formation on quasiparticle structure and the dissipation 
mechanism is strongly required. 

Recently, the carrier scattering rate in the 7r-band 
has been suggested 0, |22j to be much higher than in 
the a band. Detailed information about the quasiparti- 
cle structure and scattering rate can be extracted when 
the heat capacity, thermal conductivity and FFF resis- 
tivity data are combined. In fact heat capacity is es- 
sentially determined by the number of the quasiparti- 
cles. On the other hand, thermal conductivity is deter- 



mined by delocalized quasiparticles outside vortex cores 
because localized quasiparticles cannot carry heat. In 
contrast, energy dissipation in the FFF state is domi- 
nated by the quasiparticles localized within cores. We 
stress that all quantities, heat capacity, thermal conduc- 
tivity, and FFF resistivity in H \\ ab follow a striking sim- 
ilar trend; a crossover from low field to high field region 
occurs at H cr /H®2 —0.1 and all quantities increase up 
to nearly half of the normal state values near H cr with 
C e (H cr )/C e (H$) ~ 0.55, K(H cr )/K(H$) ~ 0.57, and 
p c f (H cr )/p c f (H^) ~ 0.45. This implies that the carrier 
scattering rate of the a band is close to that of V the band. 
This appears to be inconsistent with the suggestion by 
Ref. |2ll |23 • A more detailed microscopic calculation is 
needed to evaluate the carrier scattering rate. 

In summary, the FFF resistivity in MgB2 exhibits 
unusual field dependence, which is markedly different 
from conventional s-wave superconductors. Two dis- 
tinct regimes in which the FFF resistivity exhibits dif- 
ferent field dependence are observed. These unusual in- 
dependences indicate that two very differently sized su- 
perconducting gaps manifest in energy dissipation asso- 
ciated with the vortex dynamics. 
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